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ay: ABSTRACT 

Static stability is investigated as a contribution to numerical 
baroclinic-zone analysis, developed by Renard and Clarke, and presently 
being produced on a hemispheric basis by the U. S. Navy at the Fleet 
Numerical Weather Facility, Monterey, California (FNWF). Significant - 
level and mandatory-level data from the North American region at 002 and 
12002, 20 August 1965, are utilized in an experiment relating 850, 700, 
500, and 300-mb static stabilities and associated parameters to a FNWF 
front-location parameter proportional to horizontal shear of the thermal 
wind which is dependent on mandatory level data only. 

7 et a 

Three methods are used to compute static stability in the form of 

adn O , the most promising one involving a parabolic fit of potential 


OLnp 


temperature to the four mandatory levels. A derivative of static stability 





is used to locate baroclinic zone boundaries. Finally, the vertical gradient 
of potential temperature in combination with the slope of isentropic surfaces 
is used to determine isobaric gradients of potential temperature which 
yields an analysis of baroclinicity with resolution surpassing that of the 


present FNWF product. 
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TABLE AND ILLUSTRATIONS 


Statistical comparison of FNWF 1000-mb numerical 
fronts with NMC hand-analyzed surface fronts. 


A typical FNWF operational analysis of GG@. 


Relation of numerical fronts (GG@ ridges) with 
sea level pressure analysis. 


Vertical cross-section of average trade-wind — 
inversion phenomena over eastern Pacific Ocean /7/. 


Illustration of SS computation by Method A. 


SS analysis computed by Method A for 850-mb 
00Z 20 August 1965. 

Solid lines: SS analysis (contour interval: .05). 
Heavy solid lines: SS ridge lines. 


Same as Figure 5 except for 700-mb. 
Same as Figure 5 except for 500-mb. 
Same as Figure 5 except for 300-mb. 


Sketch of horizontal SS gradient induced by 
vertical variations. 


Schematic example of baroclinic-zone slope 
relationships. 


SS ridges computed by Methods B and C for 850-mb 
00Z 20 August 1965. 

Solid lines: SS ridge line, Method C. 

Crosses : SS ridge line, Method B. 

Numerical SS values shown above station: Method B. 


Numerical SS values shown below station: Method C. 


Same as Figure 11 except for 700-mb. 


Same as Figure 11 except for 500-mb. 
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TABLE AND ILLUSTRATIONS (cont'd) 


Same as Figure 11 except for 300-mb. 


SS analysis computed by Method C for 
850-mb 00Z 20 August 1965. 
Solid lines: SS analysis (contour interval: .05) 
Heavy solid lines: SS ridge lines. 
Crosses: FNWF o ridge lines. 
Black circles: FNWF GG8 ridge lines. 
Open circles: FNWF GG6 trough lines. 
Blocks: NMC hand-analyzed fronts 
(850-mb only). 


Same as Figure 15 except for 700-mb. 
Same as Figure 15 except for 500-mb. 
Same as Figure 15 except for 300-mb. 


SS contour chart for 00Z 20 August 1965. 
Solid lines : SS ridge lines, 850-mb. 
Dashed lines : SS ridge lines, 700-mb. 
Open circles : SS ridge lines, 500-mb. 
Crosses : SS ridge lines, 300-mb. 


Same as Figure 15 except for 12Z 
20 August 1965. 


Same as Figure 20 except for 700-mb. 
Same as Figure 20 except for 500-mb. 
Same as Figure 20 except for 300-mb. 


Same as Figure 19 except for 124 
20 August 1965. 


GSS analysis for 700-mb for 002 
20 August 1965. 
Solid lines: SS analysis (contour interval: .05) 
Heavy solid lines: SS ridge lines. 
Dashed lines: GSS ridge lines. 
Crosses: GSS trough lines. 
Black circles: FNWF GGé@ ridge lines. 
Open circles: FNWF GG6 trough lines. 
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TABLE AND ILLUSTRATIONS (cont'd) 
Page 


Computed G8 and FNWF Gé@ for 850-mb 002 83 
20 August 1965. 
Solid lines: computed G8@ contours (contour 
interval: 1°9C/100km.) 
Dashed lines: FNWF G@ contours (contour 
interval: 1°C/100km.) 
Double solid lines: computed G® ridge lines. 
Double dashed lines: FNWF G6 ridge lines. 


Same as Figure 26 except for 700-mb. 84 
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ABBREVIATIONS AND DEFINITIONS 

agebyac:eCenstants:. 

A: Degrees absolute. 

g: Force of gravity. 

FNWF: United States Navy Fleet Numerical Weather 
Facility, Monterey, California. 

GGB: -- $8, = FNWF's operational front-location 
parameter. 

GSS: A parameter similar to GG@, defined as -VSS-n,. 

G@: = ae : Magnitude of the gradient of potential 
temperature. 

Hi, The unit vector ascendant directed normal to 
isentropes on a constant pressure surface. 

NMC: National Meteorological Center, Environmental 
Sciences Services Administration, Suitland, Maryland. 

NPGS: United States Naval Postgraduate School, Monterey, 
California. 

P, p: Atmospheric pressure 

P.: Pressure on an isentropic surface. 


~) 


ean Py / Re Pressure at the first significant level above 

a mandatory level, at the first significant level below 
a mandatory level, and at the mandatory level, 
respectively. 


9S: Static stability. 


ABBREVIATIONS AND DEFINITIONS (cont'd) 
Se A ( ): A finite increment of ( ). 
16. A (AP): Error in pressure interval bounded by mandatory 
and significant levels. 
17. A(A @): Error in potential temperature change over 


interval AP. 


[fete AGB: /Computed G@ (using Equation (11)) / - /FNWF's 
processed Eny, : 

19. So. SErrorm Static stability. 

Ze ¥ : Lapse rate. 

Zale Y4 : Dry adiabatic lapse rate. 

29. go : FNWF's:static stability parameter. 

ZS 8 : Potential temperature 

24. B., OL, Om. : Potential temperature at the first significant 


level above a mandatory level, first significant level below 


a mandatory level, and at the mandatory level, respectively. 


Zo. a Slope of an isentropic surface. ' 
26. VSS: Gradient of static stability. 

Zi. V8: Gradient of potential temperature. 

28. SP : Change in potential temperature relative to a change in 


pressure at constant x, y,t. 


29. at : Change in temperature relative to a change in height 
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iL 3 Introduction 

Considerable developmental work has been performed at the U. S. 
Navy Fleet Numerical Weather Facility (FNWF) and the U. S. Naval Post- 
graduate School (NPGS) on the problem of locating frontal zones by 
numerical methods. W ~o" 9 / The most promising and currently operational 
approach utilizes the second derivative of potential temperature (- a2 2 ) 
as a front-location parameter, hereafter referred to as GG@. Axes of relative 
minima and maxima in the GG@ fields which bound the baroclinic zones 
are associated with relative maxima in the horizontal shear of the thermal 
wind. The axes of relative maxima represent the warm-air boundary of 
baroclinic zones and, as such, are a numerical form of conventional 
fronts. These numerical fronts have a close resemblance to hand-analyzed 
fronts, as may be seen in Figure 1 which shows a typical FNWF operational 
analysis of GG@, and Figure 2 which indicates the relation of the numerical 
front (i.e. , GG® ridge) to the sea-level pressure analysis. 

However, recent statistical evaluations indicate a need for-improving 
the numerical product in certain areas, such as in eastern ocean and sub- 
tropical regions and frontal occlusions. The author of this study compared 

1throughout the subsequent discussion the following terminology 
is used interchangeably: (a) GG@ fronts, numerical fronts, baroclinic 
zone boundaries; (b) GG8@ ridge, numerical front, warm boundary of 
baroclinic zone; (c) GG®@ trough, cold boundary of baroclinic zone; 
(d) numerical frontal zones, baroclinic or hyperbaroclinic zones; 


(e) hand-analyzed fronts, manually analyzed fronts, conventional 
fronts, NMC fronts. 


ll 


1000-mb fronts, specified by GG@, with hand-analyzed surface fronts of 
the National Meteorological Center, (NMC) for 16 observation times in the 
late summer of 1965 ew. It was found that 72% of the GG@ fronts could 
be associated with NMC fronts. The average distance separating the two 
fronts was about 90 nautical miles, slightly greater in sparse data areas. 
These statistics were based on a sample of over 250,000 miles of fronts. 
Table I recapitulates this and other findings. 

In its present Operational form, 8 is determined from mandatory level 
data only, however 89 or its derivative (G@ and GG8@) may be computed along 
any level of meteorological interest (e.g., 1000-mb, earth's surface, 5,000 
feet, mean level of 1000/700-mb layer.) It is reasonable to expect that 
consideration of the detailed temperature structure (i.e., from significant 
as well as mandatory levels) in the vertical will contribute to the accuracy 
of locating baroclinic zones at the aforementioned levels. 

This study introduces detailed vertical resolution through a static 
stability parameter computed from significant-level as well as mandatory- 
level data. The quantity was evaluated, first as an entity unto itself, and 
later, combined with other quantities to determine frontal locations and 
intensities. Specifically, the computed static stability field was considered: 

(1) Relative to FNWF's operational model of GG@, as well as to their 
present static stability evaluation /5 : ee (See section 4) 

(2) As a contribution toward a more detailed GG@ analysis through 
the introduction of detailed vertical temperature structure. (see section 6) 
on Selection of a Suitable Static Stability Parameter. 


Static stability refers to the resistive force offered to motions of 
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various types by the vertical density structure of the atmosphere in the 
gravitational field and, specifically, for the atmosphere in hydrostatic 
equilibrium. Here static stability is proportional to Ys -¥ . When Py , 
an individual air parcel displaced dry adiabatically in a vertical direction 
will tend to return to its original position, and the air column is said to be 
statically stable. 

Frontal layers are zones of relatively high values of static stability, 
since within such layers temperature usually deviates from its normal 
decrease with increasing height to yield isothermal or inversion conditions. 

Thus, it appears meaningful to compute static stability for geometric 
or isobaric levels of interest over an area. Relatively high values of static 
stability ought to be located within frontal zones and, thus, the static 
stability should have a high positive correlation with baroclinicity; GG®@ 
ridges and troughs should bound zones of high static stability. Since a 
large value of stability does not imply the existence of a baroclinic zone, 
and vice-versa, it is expected that some stability ridges (that is, axes of 
relative maxima) and numerical fronts may exist independently. Figure 3 
shows a low tropospheric vertical cross-section of the eastern Pacific Wa7/. 
This cross-section illustrates high values of static stability and baro- 
clinicity in the east Pacific trade-wind inversion east of 135W, while also 
indicating an area of poor correlation between baroclinicity and stability in 
the inversion west of 135W. 

Gates has made an analysis of advantages and disadvantages of 
several forms of static stability /8/. Selection of an optimum stability 


parameter for this project was based on several self-imposed requirements: 
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(1) static stability will be computed on constant pressure surfaces’; 
(2) the measure should have equal, or nearly equal, values at different 
levels when ¥ is the same fraction of al : (3) static stability should 
be zero when & = Xd - (4) input parameters should be readily available 
from radiosonde data; and, (5) the computation must readily lend itself to 
both manual and numerical analysis, the latter currently existent or planned 
for the future. Requirements 1, 4, and 5 are self explanatory. Requirements 
2 and 3 are necessitated by a desire to have comparative fields of static 
Stability at the levels investigated. 

A relatively simple, non-dimensional measure possessing the desirable 
properties is: 

$$ =-} = (1) 

This measure meets the above requirements very well, moreover, it-is 
directly proportional to Ye -% . Gates found that SS computed by this 
method showed a slight decrease with. height up to 300-mb, followed by an 
increase, but the range of this variation is small and will not be of any 
consequence here. Equation (1) also portrays the stratosphere realistically, 
and shows the relatively unstable upper troposphere in both summer.and 
winter. 
3. Area and Time of Experiment. 

The area and time selected for computing SS were all of the United 
States, Canada, and Alaska at 00Z and’ 12g; 20 August 1965. The reported 

2rhe present FNWF operational frontal model utilizes mean 9 and 


its derivatives for the layer 1000/700-mb as the basic frontal parameter. 
Derivatives are computed at constant pressure. (See Figures 1, 2) 
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radiosonde temperature data were then utilized in three different methods 
to compute values of 5S, for each of approximately 100 reporting stations, 
af's50, 700 , SOO*"and*"S00=nip” 

20 August 1965 was a well documented day, with extensive data 
readily available from listings provided by FNWF. The date was marked 
by the presence of a strong cold front and large stability over the eastern 
and southern United States. As Gates pointed out, wintertime situations 
give greater stability contrasts (up to 50%), especially in the lower tropo- 
sphere, and better delineate frontal zones, but 20 August proved to be quite 
satisfactory for a pilot project. 

4. Computation and Areal Analysis of Static Stability. 

I. Methods for Computing Static Stability. 

a. Method A. This method consists of computing the stability 
for the two significant [aPePSadlacerit to a mandatory level and plotting 
the larger of these two values on the appropriate constant pressure chart 
(850, 700, 500, and 300-mbs.) This approach has the desirable effect of 
sharply defining the stable baroclinic zone boundaries by assigning to the 
mandatory level the highest possible SS value from one of the two adjacent 
significant layers. Figure 4 shows a portion of a typical sounding for which 
SS might be computed using Method A. Since layer AB is more stable than 
layer BC, SS is computed for AB and-this value is assigned to 500-mb, the 
mandatory level. Here, the 500-mb level represents the warm-air boundary 


3A significant layer is bounded by two significant levels, one or 
both of which may be mandatory levels as well. 


IS 


of a stable and, presumably, baroclinic zone. Thus a maximum stability 
value found in a layer within the stable zone is placed at its boundary. 

Stability within the significant layer was computed according 
to the definition in Equation (1). P and @ are the median values for the 
layer, while the layer value of os is approximated by BP of the layer. 
In this method, a significant layer was defined as being bounded by a 
mandatory level and its nearest significant level, regardless of the 
layer's depth. In another method to be discussed later, layers of less 
10-mb thickness were ignored. 

A Fortran 60 SOG was written to compute stability to four 
significant figures by this method. The values were rounded to three 
significant figures and plotted and contoured on FNWF's 1:30 million polar 
Steoreographic projection. Resulting analyses appear as Figures 5-8, with 
comments on the technique appearing in Appendix I. 

Figure 6 (see station 206) displays the merit of this approach, 
since large gradients of stability do allow the stable zone boundaries to 
be easily located in this case. 

But Method A was not considered satisfactory for several reasons. 
First, it must be remembered that the objective is to compute stability at 
the mandatory levels. This system arbitrarily assigns a layer value as the 
mandatory-level value. The actual value at the mandatory level approaches 
this arbitrary value as a limit, but for many stations this method gave the 


4This and all subsequent programs were run on the NPGS's 
Control Data Corporation 1604 digital computer. 
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stability for a small layer near the mandatory level, which is not nec- 
essarily representative of stability_at the mandatory level. 

Secondly, this method introduced vertical variation into hori- 
zontal gradients of SS. Figure 9A shows a hypothetical case in which SS 
is higher in the significant layer above station A than in the significant 
layer below station B. The resulting field of SS tends to give excessive, 
unrealistic gradients of SS within the baroclinic zone, as well as suggesting 
improper boundaries of the zone. This problem is magnified for stable zones 
which are thin (less than 20-mb in depth) and undergo marked variations in 
depth over small horizontal distances. 

Thirdly, it will be noted from a discussion of Methods B and C, 
that the latter give a more meteorologically satisfying field of SS relative 
to GG@ baroclinic zones and the hand-analyzed fronts of NMC. 

Consequently this method did not look as promising as Methods 
B and C, and was not further developed. | 

b. Method B also utilized Equation (1) for each of the four 
mandatory levels of interest by employing a mathematical best fit of 9 at 
the mandatory level, using transmitted radiosonde data. This data consists 
of temperature and pressure recorded at each significant level (level where 
a significant change in temperature lapse rate occurs), as well as for all 
mandatory levels. Thus the sounding, as plotted from transmitted data, 
is discontinuous at each of these levels. However, analysis of radiosonde 
data, aSsuming every temperature contact is a significant level, indicates 


the transmitted discontinuity points are largely manufactured by coding and 
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processing procedures. Thus, fitting a continuous curve to the transmitted 
sounding for the purpose of obtaining oy is reasonable and justified /2/. 


This is accomplished by numerically solving the following system of 


simultaneous linear equations for each sounding: 


naa bpu + CDae 
OmL= at bDmni tC pu (2) 
Ov = ORE DP +CDi~ 


A Fortran 60 program using a Jordan 2 subroutine was written to 
solve this system for the constants a, b, andc, Since all six pressures 
and temperatures (hence, 9) are known from the sounding. 

At the mandatory level, the vertical gradient of potential tem- 


perature, 22 becomes 


ap 


29m = b+ 2C- pm . 


Thus, Equation (1) was numerically evaluated for 
Saat ee a ee (Dt CP) (4) 


It can be shown (Appendix III) that the finite difference form 


ale) 
of 5 Prot 


is: 


ap Lape AOL + RES x Ou | (5) 





This expression was used to manually check computerized 


results of Equation (4). One such computation appears in Appendix II. 
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For FNWF numerical operational purposes, Equation (5) would undoubtedly 
prove to be a more desirable system than Equation (3). 

Computed static stabilities were plotted and contoured ata .05 
interval, similar to values resulting from Method A. 

Two problems appeared at this point: (a) the existence of some 
negative values of SS on each chart (i.e., absolutely unstable conditions) 
which verified to be real and representative of low level Summer conditions, 
and,thus, were accepted; and (b) a few questionable values which caused 
irregularities in otherwise smooth contour patterns. Re-examination of 
sounding data revealed that a majority of these questionable values of SS 
resulted when a Significant level was located within 10 mb of a mandatory 
level. 

When thin layers of a few millibars are considered, normal 
errors in pressure and potential temperature can vary the value of SS at 
the mandatory level by a large percentage. This problem is discussed more 
fully under Method C. 

Furthermore, layers of a few millibars' thickness may be con- 
sidered mesoscale in character and not sufficiently deep to define synoptic 
scale frontal layers, or, alternately, as nearly horizontal, thin, stable 
layers whose relation to baroclinic zones is unexpectedly poor. For 
example, near 700-mb, a 10-mb layer (about 400 feet) representing a 
narrow baroclinic zone 100 miles wide, yields a 1/1500 slope. Since 
baroclinic zones as depicted by GG@ average almost 300 miles, a thin 


layer (10 mb or less) is, for all practical purposes, horizontal. (Figure 10) 


IES, 


By ignoring layers of 10 mb or less, this system would be throwing out 
only the nearly horizontal stable layers. 

A third method has now been suggested, one which does not 
consider layers of less than 10-mb thickness adjacent to a mandatory level. 
This is Method C, discussed below. Since Method A altered 58 out of a 
possible 440 values computed in Method B, taking all four levels into 
account, the location of SS ridge lines differed for each method, but contour 
patterns were basically the same. For this reason, the contour analysis by 
Method B was not included in this paper. Instead, Figures 11-14 compare 
the locations of the SS ridge lines produced by Methods B and C. These 
figures are discussed in Appendix II. On each chart, SS values are shown 
for all stations where SS varied depending upon the method used. The 
number plotted above the station is the value computed by Method B; the 
number below the station is Method C's result. It is evident that the choice 
of methods heavily influences the location of ridge lines, and the choice of 
a minimum significant-layer thickness can be a critical factor in locating 
SS zones. 

c. Method C. This method is considered as best representing 
SS at the four levels of interest. It is identical to Method B, with one 
exception. Layers 10 mb and less in thickness are considered unrepresent- 
ative. In cases where a significant level exists within 10 mb of a mandatory 
level, the next lower or higher significant level is substituted. For example, 
assume the levels of interest are 886, 858, 850, and 795 mb. For 850-mb 


analysis, 858 mb would be ignored since it is within 10 mb. Instead of 
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858 mb, 886 mb and its associated temperature are substituted for a 
ander " 
u 

Thus , Method B is a very precise method, taking every fluc- 
tuation of temperature in the layers surrounding the mandatory levels into 
account. Method C includes a small "Smoothing factor", which has the 
effect of filtering out some of the "noise" created by very thin (i.e., non- 
synoptic scale) layers. Furthermore, the minute temperature changes in 
thin layers dictates use of layers thicker than 10 mb, primarily because of 
the accuracy demanded by this method. 

When compared with FNWF numerical fronts for the same time 
and level, the more satisfying appearance of the SS ridges using Method C 
raises the question as to the best smoothing factor to be employed. Per- 
haps raising the minimum thickness of the significant layer to 20 mb may 
even be justified. An error analysis of SS brings out this point. Consider- 
ing an error in static stability (ASS) at a mandatory level as a function of 
errors in Significant layer depth Vix (A P) /, and temperature change in the 


layer Van (A WA Equation (1) may be differentiated to obtain: 
- -| 
ASSme = -[-(A8)(A2)* A(AP)+AP A(A@y] Lee 6) 


Probable errors in pressure and temperature are assumed to be 1 mb and 
19°C, respectively, in the low troposphere Par/ So as not to inflate the 
proposed error, it may be assumed the probable error occurs at only one 
boundary of the significant layer, errors in 8 and T are the same magnitude, 


and transmission errors are excluded. Further, assuming P = 700-mb, 
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9 = 310°, and AQ@= .2°C for each AP interval of 2 mb, the following 


results: 
A P(mb) A6(°C) ASS 
2 2 104% 
10 li 21% 
20 2 11% 


Thus, dramatic changes in accuracy are achieved by proper 
selection of the smoothing factor, while still retaining the desired vertical 
resolution. Further experimentation is certainly warranted before final 
selection of the optimum minimal thickness of significant layers. 

d. Analysis Techniques Employed with Method C Data. 

Data obtained by Method C were plotted on 1:30 million 
polar stereographic charts and analyzed for 00Z and 12Z 20 August at all 
four levels. (Figures 15-24) 

Any attempt to contour an area as large as the United States, 
Canada, and Alaska on the basis of about 100 observations is open to 
criticism, since any one of several interpretations may be used without 
violating data. This is especially true in Canada, where extensive sparse 
data areas exist. The proper orientation of the isolines of SS is intrinsic 
to proper location of the ridge line of SS. Since the ridge line suggests the 
location of the frontal zone and/or its boundaries, an assist or first guess 
was sought for the SS analysis. 

Since the ultimate goal of this project is to guide, or even 
replace, the present FNWF system of frontal analysis, it was considered 
valid to use the GG@ field as an aid in contouring SS by orienting the SS 


lines parallel to the numerically computed baroclinic zone boundaries 
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wherever data allowed. It turned out that forcing was not necessary, 
since SS isolines did indeed parallel baroclinic boundaries on most 
occasions, especially in areas of strong baroclinicity. The result looks 
very promising, since the patterns obtained show a striking resemblance 
to GG@ ridge lines, notably at 850 and 700 mb. (Figures 15, 16, 20, and 
21) 

II. Discussion and Evaluation of SS Charts Produced 

At this juncture it was felt that Method C, considering baroclinic 
zones presented SS in its most palatable form as representative of frontal 
phenomena. The only obvious shortcoming is its subjectivity, being a 
manually produced analysis product. 

Appendix III contains a detailed analysis of eight charts, one 
for each level of interest for both 00Z and 1200Z 20 August 1965, the charts 
appearing as Figures 15-18 and 20-23. Speaking in general terms, it can 
be observed that: 

a. Most ridge lines of SS are in close proximity to NMC fronts 
and to ridge lines of GG@; exceptions are manifest at higher levels (500 and 
300-mb), and in mountainous areas. 

b. There are more miles of SS ridge lines than GG@ fronts, and 
more GG@ fronts than NMC fronts, the latter fact indicated in vaio. and yids / | 

c. SS ridge lines have obvious vertical consistency between the 
850 and 700-mb levels, especially over the dense data, non-mountainous 
region of North America, where a hand-analyzed front also exists. Here, 


the 700-mb front trails its 850-mb counterpart by about 100 miles. This 
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seems consistent with the normal slope of a low tropospheric front. In 

the same area, the 500-mb SS ridge trails the 700-mb position by as much 
as 300 miles (Great Lakes area, 1200Z) and it is at this level that vertical 
continuity breaks down. Not all fronts reach the 500-mb level, and of those 
that do, many show a sharply decreased slope. Such an observation is in 
agreement with findings by other investigators VOT Vertical consistency 
of SS ridges in mountainous regions, sparse data areas, and the peripheral 
map area is fair to poor. 

d. The ridge lines at 300-mb are questionable representatives 
of low troposphere fronts for the same reasons outlined above. There is the 
additional problem posed by the tropopause. In computing SS by Method C, 
data up to 250-mb were often used. In winter, the low-level tropopause 
gives high stabilities over some mid-latitude stations, and thus gives a 
false indication of a front. Several soundings in Canada (stations 722, 
909, 913, 915, 917, and several others) indicated a tropopause existing 
between 250-éand.300-mb 

e. The ridge lines of SS showed good time continuity in the area 
of good vertical consistency. 

f. In areas of Sparse data, subjectivity enters the analysis. In 
areas where three or four reports determine the contour pattern over an area 
of thousands of square miles, baroclinic zones, history, and vertical con- 
sistency are the only aids etapa. These are important considerations 
when viewing application on a hemispheric scale, utilizing computerized 


solutions. 
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III. Comparison of SS with FNWF's Static Stability Analysis (0 ). 

a. Observations: 

FNWF's data processing program produces mean static stability 
for five layers (0 ) / S86/. Such a computation, of course, was not meant 
for aid in frontal analysis. However, it is pertinent to this investigation to 
compare the merits of 5S computed for a level, using the maximum vertical 
resolution, with 0, computed for a layer using mandatory level information 
only. A detailed comparison of the two analyses appears in the comment 
section of Appendix IV. The discussion here is limited to the following 
general observations. 

The use of mandatory level data only gives a gross picture of 
stability, and says little about stability at a particular level. A frontal 
layer could be embedded within the 775/600-mb layer, and mandatory level 
data might never detect it. On the other hand, 5S using significant as well 
mandatory level data is more representative, if for no other reason than the 
increase in data and the vertical resolution. 

5S has the added advantage of being computed for four specific 
constant pressure surfaces in the atmosphere. With the wide use of 850; , 
700, , and 500-mb charts, stability computed at a level is more meaningful 
than stability within a gross layer. 

b. Conclusions: 

When compared with FNWF's © analysis, 55 is a more 


meaningful representation of frontal zones because of: 
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(1) the use of actual data from radiosonde reports vice 
processed data as used by FNWF. 

(2) the high vertical resolution resulting from the use of 
significant vice mandatory level data. 

(3) the improved accuracy derived from using a parabolic fit 
of potential temperature at the mandatory level with no assumption about 
the vertical gradient of static stability. 

IV. Recommendations on Ss. 
Based on the aforementioned findings, the following recommenda- 

tions are made relative to further testing and application to frontal analysis: 

a. Limit the mandatory levels to the 850,, 700, , and 500-mb 
surfaces for low tropospheric frontal phenomena. 

b. Compute SS as indicated by Method C. 

c. Statistically compare SS ridges with computer and hand- 
drawn fronts and baroclinic zones produced by NMC and FNWF, respectively. 

d. Consider programming computation of SS on a hemispheric 
basis, say in the framework of the data processing scheme at FNWF. It is 
realized that this requires, apriori, a Suitable data-handling program in 
order to marry the mandatory and significant level radiosonde data. Pre- 
sently, the significant level data are copied on tape but not read out or 
utilized in any way in FNWF's automatic data processing program. 

e. Subsequent tod, ee the SS field considering it or its 


derivatives as an adjunct to the GG@ analysis for locating baroclinic zones. 
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f. Investigate the feasibility of computing SS for the 
surface, or 1000 mb, by assigning the SS computed for the lowest 
significant layer (or level) as the value at the surface. 
om The Gradient of Stability: GSS. 

I. Definition of Parameter. 

As indicated in the Introduction, FNWF's present operational 
product for frontal analysis is the GG@ field. This field, as defined by 
Renard and Clarke, is the directional derivative of the gradient of potential 
temperature along its gradient Ve In mathematical terms: 

o* oO 


GGO = -V/VO/l WB -= -V/Vel- eS = (7) 


a) 


Referring to SS and recalling that it is already a first derivative, 


a similar expression involving stability may be written: 





GGS= -VSS YS PT GGe Te (8) 


The GG@® and GSS fields are entirely similar in gross features,as 
may be seen by Figure 25. 
II. Discussion of Procedure. 
GSS was calculated manually for the 700-mb level for 002 
20 August 1965. The mechanics of this calculation were as follows: 
a. Using FNWF's hemispheric printout of potential temperature 
at 700-mb, re was determined at each grid point over the area of interest. 
b. For each grid point, the gradient of stability was deter- 
mined using Ao> , where An, was taken as 120 nautical miles, with 
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the grid point at the centrum. This distance is rather arbitrary, being 


Ld, 


selected mostly on the basis of allowing for a well defined field of GSS. 
In other words, the larger the distance Ang: the flatter the field of GSS. 

c. GSS was plotted on a 1:30 million polar stereographic 
projection and then analyzed. This chart appears as Figure 25. The warm 
air boundary of a baroclinic zone is associated with the ridge or axis ofa 
relative maximum in GSS, while the cold air boundary is representative of 
the trough or axis of a minimum is GSS. 

lL Genelusions . 

a. GSS does appear to offer some advantage over the SS field 
for frontal information. SS defines a ridge only, which is associated with 
the centrum of the zone, while GSS defines a ridge and a trough which out- 
line the boundaries of the baroclinic zone. In comparing these boundaries 
with the baroclinic zone limits derived from GG@, the GSS boundaries define 
a zone which averages about 100 miles narrower than GG8@, except over the 
west central U. S. (See Figure 25) There are at least two reasons for the 
suggested narrower baroclinic zones: 

(1) greater vertical resolution in computed GSS vs. GG@ 
from FN WE, or 

(2) use of 120 miles to compute GSS instead of several grid 
lengths, as used by FNWF. 

b. The manual computation of GSS is a time consuming process, 
requiring about three hours to Ereanee the 700-mb chart for the North 
American area. However, the computation of GSS is adaptable to the 
computer, the mechanics of the operation being almost identical to the 


computation of GG@. 
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IV. Recommendations. 

a. If GSS is investigated at some later date, the procedures 
outlined in this section are recommended, but it is suggested that the 
entire operation be programmed with GSS as the final product. 

b. If the above recommendation is carried out, statistically 
compare GSS with GG@ to determine if the former is consistently better 
in outlining the baroclinic zones. 

c. The following section concerns an alternate approach to 
obtaining baroclinic zone boundaries. Since SS and GSS experiments were 
but pilot projects, it is recommended that a more exhaustive study be made 
relating the GSS scheme to the approach in the next section. 

Ge. G8 Computed Using High Vertical Resolution. 
I. Definition of Parameter. 

As discussed in Section 3, GG8@ in its present operational form 
utilizes mandatory level data only for frontal analysis. 

Consider the following: 

If O=¢46n. (P,Ne) 

dO= 3%, 6No+ 5% Sp (9) 


Next, set d@-O for purposes of computing the slope of an 


isentropic surface, here, in the direction of maximum slope. Thus, 

















2 OU 20 
SNeale a * (10) 
and, Se P| (11) 
where, GO= ——" 2) 
Binelly weeci—e-GiGe) = ———> 
t are 
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II. Discussion of Procedure. 

Equation (11) can be solved from available data since 2°. 
may be computed using Method C, and SF. can be determined 
utilizing FNWF's hemispheric printout of eat Thus, a G8 results which 
incorporates the high vertical resolution inherent in the use of significant 
level data through consideration of 3p . This new G@ may then be com- 
pared with FNWF's hemispheric printout of G@ to determine whether high 
vertical resolution used in this fashion enhances the operational frontal 
preduct™ 

The mechanics of determining bP | involved obtaining hemi- 
spheric printouts of P for each 2.5 A of potential temperature from 280A 
to 330A, in order to compute - 9 at 850-' and 700-mb. for 00Z 20 August 
1965. | Wale was computed at each station, using a centered finite 
difference scheme applied to two mesh lengths of the FNWF grid. Inter- 
polation of ace for the exact 8 surface was made, as appropriate. 
Finally, the products of ane and ys were Obtained for each station 
at 850; and 700-mb. The resulting values of G@ (hereafter referred to 
as computed GQ) were transferred to 1:30 million charts and compared with 
FNWF's G8 values. 

III. Results. 

The comparison between computed and FNWF G6 values appears 
in tabular form as Appendix V. Note that computed G@ has been uSed as the 
basis, Since its values are correct to three significant figures, and its use 
allows the selection of 177 specific reference points (considering both charts) 


for comparison of the values. 
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Each chart was arbitrarily separated into two parts, for reasons 
which will soon become evident. One part was the baroclinic zone, which 
was defined here as the area enclosed by the 1.0 C/100 km. iselime.> The 
remaining area was Simply considered outside the baroclinic zone. Figures 
26 and 27 show a comparison of the baroclinic zones from computed and 
PNWF G6 values. 

The table in Appendix V reveals that computed G® gives relatively 
higher values in baroclinic zones and lower values outside the zones. This 
fact is obvious when one notes that of 51 reference points within the baro- 
clinic zone, G@ is less than ~ .2°C/100km at two points only; while of 126 
reference points outside the baroclinic zone, only three values of AG®@ 
are above +.2°C/100km. 

Thus, as shown by the medians obtained for each area, computed 
G® gives a wider range of values overall, implying a sharper change in 
gradient (i.e. , GG@) across the boundaries of the zones. 

The 16 observations at 850-mb in which the FNWF G® exceeded 
computed G@ by more than 1 .0°C/100km fell mostly in the southwest desert 
area of the United States, as may be inferred from Figure 26. As stated in 
Section 1, this area poses a problem for numerical frontal analysis by GG@ 
in that fronts are indicated where none exist in manual analyses, especially 
in the Summer. Computed G@ has markedly low values in this area. Thus, 
either the vertical resolution of aye or the use of a three-point vice a 

attics onlveiki charts ranged from nearly 0 to 3.0°C/100 km. 


with two isolated higher values on computed G@: 5. 15°C/100km at station 
764 and 3.50°C/100 km at station 518. 
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five-point derivative in computing xa has, at least for this map time, 
resolved an area of difficulty. 

A final result concerning locations of relative maxima and 
minima in the field show that of 23 relative maxima on the computed G@ 
field, there were 17 maxima on FNWF's printout which fell within one mesh 
length. Of 29 relative minima, FNWF's printout indicated 23 minima within 
one mesh length. These facts, combined with a cursory glance at Figures 
26 and 27, show a Similarity of contour patterns, except over the southwest 
Ur Swat 850-mbr" 

IV. Conclusions. 

It is difficult to make any sweeping statements about G@ 
computed by Equation (11) on the basis of one map time, but the following 
comments certainly appear justified: 

a. There is a reSemblance in the shape and location of 
baroclinic zones defined by the two methods, especially the warm air 
boundary. 

b. Computed G@ defines a narrower baroclinic zone with 
higher values inside and lower values outside the zone. 

c. Computed G@ seems to eliminate the troublesome area 
Over the mountainous and desert southwest at 850- and 700-mb. 

Thus, computed G®@ appears to show considerable promise, 
especially in light of b. and c. The question reduces to an analysis of 
Equation (11) and a decision as to whether the vertical resolution obtained 


by the use of significant level data has narrowed and sharpened the 
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baroclinic zone, or whether the use of a three-point vice a five-point deri- 
vative is responsible. Comparisons made by Associate Professor R. J. 
Renard of the U.S. Naval Postgraduate School and Leo C. Clarke of 
FNWF, Monterey, California, revealed that the five-point method actually 
gives a stronger gradient than the three-point. Thus, it certainly appears 
that Bie plays a dominant role in statement b. 

As a test, three stations at 700-mb were arbitrarily selected 
and a comparison made between the three-point (two mesh lengths) 
derivative used in this experiment and the five-point technique presently 
used operationally by FNWF. The five-point derivative is described by the 
quartic interpolation polynomial 

OP av i= O= aaqg € (Plia + BLOP)i+ - (Pin - (P)iv23 (13) 
Where d= mesh length = 38lkm. at 60" latitude, i= 0 is the station 
location, i+ 1 and i+ 2 are locations one and two mesh lengths from the 
station, respectively, in a direction normal to isobars on a constant 8 
surface. The table below suggests that the gradients computed by the two 
methods are similar, but certainly not the same. 


&p/ § n,(mb/100km) 


Station three -points five-points 
7a 32.8 AQ v3 
as 40.6 4 Seed 
oe Lied 14.6 


V. Recommendations. 


Continue experimentation on computed Gé@ by: 
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a. programming the computer to produce 26 from Le 
charts and multiplying the result by 25 from Method C; 

b. operate on computed G89 to get GGO (this is already an 
operational procedure at FNWF) with high vertical resolution; 

c. make a detailed statistical comparison of the two GG@ 
parameters to determine which product best delineates the baroclinic zone 
and its associated frontal phenomena. 

7. General Conclusions and Recommendations. 

In view of the stated conclusions in each section, only a brief overall 
review will be given here. 

This entire study utilizes data from two map times and attempts to 
show the value of using significant as well as mandatory level data as a 
contribution toward numerical baroclinic zone and frontal analysis. Three 
different parameters for frontal analysis have been presented: SS. GSS, and 
a computed G8. It is not the purpose of this investigation to select from 
these three a single proposal which is considered best. All three have 
strong and weak points, anda statistical study covering several periods of 
time in at least two seasons would be required to make an intelligent selec- 
tion. 

It does appear, however, that the detailed thermal structure of the 
atmosphere in the vicinity of the mandatory levels can be better visualized 
by use of significant level data, and that the computation of oa using a 


r9) 


parabolic fit is a most promising approach. 
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TABLE 1 
EVALUATION ANALYSIS 


FNWF (NUMERICAL) 1000-MB vs NMC (MANUAL) SFC FRONTS 


Dates: 16 obsn. times in period 23 Aug - 6 Sept 65 
FNWF FRONTS: 259,000 n.mi. 
NMC FRONTS: 284,000 n.mi. 
% FNWF FRONTS assoc. with NMC FRONTS: 72% 


% NMC FRONTS assoc. with FNWF FRONTS: 66% 


Frontal Average Difference (lat) 
Sector Distance (FNWF vs NMC) 

(n.mi.) Summer Winter 
North America 49 ,250 1.6 1.4 
Atlantic 46,900 1.4 1.8 
Pacific 46,300 Aol Bes 
Europe-North Africa 26,000 1.8 --- 
Asia 18 ,000 ee, =<-- 


een VP /) 
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APPENDIX I 
SS Analysis by Method A 
Technique and Criteria. 

SS computations by Method A were plotted on 1:30 million polar 
stereographic paper and the resulting contour analysis appears as Figures 
5-8. Values of SS vary from slightly negative to +.6, with most values 
falling in the +.05 to +.3 range. Contour interval has been arbitrarily set 
at .05 for values from 0 to .2, and at .1 thereafter. 

Contours of SS are generally elongated and elliptical, somewhat 
resembling a jet stream analysis. The elongation generally appears 
parallel to ridge lines. 

Ridge lines are drawn only for relative maximum values greater than 
.1 and then only if more than one station shows such a magnitude. Ridges 
are also omitted for areas of weak gradient (although values of SS may 
exceed .1) as well as when considerable subjectivity is required in locating 
the ridge. Although this discussion concerns Method A, the principles 
apply to every analysis attempted in this experiment. 
850-mb 00Z 20 August 1965. (Figure 5) 

It will be recalled from Section 4 that Method A has the advantage 
of producing numerical values of SS which are somewhat higher than those 
computed by the other methods discussed herein. These higher values have 
produced a ridge over the eastern U. 5. which is well defined by high SS at 
Stations 712, 528, 520, 429, and 403. This ridge has not been connected 


to the ridge over the central U. S. because of low SS at station 327, buta 
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slight relaxing of criteria could have allowed such a marriage. The other 
three ridges on this chart are a function of a very few station reports, and 
exact placement was difficult. 

The western U.S. has few contours because many of these stations 
are at pressures greater than 850-mb. The zero contours in the U.S. 
southwest desert region bounds areas of absolutely unstable lapse rates. 

The overall pattern has some resemblance to the frontal patterns 
achieved by other methods (see Figure 15), but it is clear that correlation 
is too poor for a comparative analysis to be of any value. 
700-mb 00Z 20 August 1965. (Figure 6) 

This chart is characterized by several high values of SS which dictate 
the location of ridge lines. Most of these values are computed from thin 
layers adjacent to the 700-mb level, and the overall pattern shows very 
poor correlation with patterns produced by later methods and operational 
products. (see Figure 16) The most striking example appears at station 206. 
Other computations modify this 1.86 reading to a value of about .1, a 
remarkable change. The same criticism can be made of stations 722, 907, 
endeS ban 

With the exception of the ridge over the eastern U. S., vertical 
continuity with the 850- and 500-mb levels is poor. 
500-mb 00Z 20 August 1965. (Figure 7) 

The 500-mb chart, when compared to GG@ and later methods in this 
study (Figure 17), has the highest correlation of any chart of the series, 


showing a long ridge of SS from SE Canada, across the centralU.S., 
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then north into central and northern Canada. High SS values at stations 
520, 712, 747, and 913 dictate the position of the ridge, just as they do 
in Figure 15. The ridge over the eastern seaboard is well defined, but 
its position is not verified by any other method. A large portion of the 
western U.S. exhibits moderately high values of SS, but the area is so 
broad and the gradient so flat that it is not possible to draw a realistic 
ridge line. 

300-mb 00Z 20 August 1965. (Figure 8) 

The pattern at 300-mb is random and irregular. The only ridge 
showing vertical continuity with lower levels is the ridge over eastern 
Canada and the east-central U. S. As indicated in the text, there is 
evidence to support the claim that the 300-mb analysis may not be an aid 
to lower troposphere frontal analysis. Its value as a front-locator is, 


therefore, questionable. 
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APPENDIX II 
Comparison of SS Ridge Lines, 
Method B vs Method C 
Technique and Criteria. 

The general discussion of Methods B and C points out that the 
techniques are the same except that in Method B all significant layers were 
considered. 

Figures 11-14 point out the differences in location of SS ridges caused 
by differences in SS values at individual stations where thin layers (=10mb) 
were present. Each analysis was accomplished without referring to charts 
produced by the alternate method, and therefore are as objective as possible. 

Criteria for ridge location are the same as those discussed in Appendix 
I for Method A. 
850-mb 00Z 20 August 1965. (Figure 11) 

This chart illustrates how changes in two values (stations 327 and 
451) can adjust the location of SS ridge lines by as much as 200 miles. 

The lower value of .000 calculated by Method B forces the ridge well north 
of station 327, and away from the NMC front through the area (Figure 15). 
The other ridges are identical, since no other thin layers were encountered. 
700-mb 00Z 20 August 1965. (Figure 12) 

Numerous thin layers at 700-mb force a Method B pattern which differs 
considerably from Method C and GG@ ridges as shown by Figure 16. As can 
be seen in that figure, Method C ridges have very high correlation with the 


FNWF numerical fronts, but Method B produces a wavy, random pattern 
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which does not correlate well with GG@. It appears that the elimination 
of thin layers (Method C) has smoothed and adjusted Method B into a 
meaningful SS ridge analysis. Certainly this 700-mb chart seems to 
justify the arguments presented in section 4. 

500-mb 00Z 20 August 1965. (Figure 13) 

The 500-mb analysis shows, once again, how changes in SS values 
tend to alter the SS ridge lines. Since verification of both Methods B and 
C with FNWF's numerical fronts is fair to poor on this chart (see also 
Figure 17), comments will be limited to pointing out that Method C again 
seems to smooth out the ridge pattern, yielding a more meteorologically 
satisfying pattern. 
300-mb 00Z 20 August 1965. (Figure 14) 

Method C simply modifies a "no correlation" pattern to another 
"no correlation" pattern over central Canada. Only five values differ, and 
the pattern by both methods is subjective due to sparse data and a wide 


variation of values. 
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Derivation of Equation (5) 


1s G=at+bpt+cp- 
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As an example of this equality, consider station 520 


Pittsburgh, Pennsylvania, at 850-mb. 


PRESSURE P. TEMPERATURE 
883 298.6 
850 pOO.1 
771 301.8 


Computer result, using parabolic fit, SS = .1082 





By tanite diftiercncing : 


ao — 1 fA sa0) Ap 10, | ‘ 
OP APL APY fopgeel 
ae. i eS eee: in7w= OnE 
West Resa 1 ZS 
SS =p20 = 850 x .0383 = .1083 
83D SO aul 


Error in Lourth@place. 
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APPENDIX III 
SS Analysis by Method C 
Comparison of 5S with GG@8@ and NMC Fronts 
Technique and Criteria. 

SS computations by Method C were plotted and analyzed applying the 
principles discussed in SS Analysis by Method A. (Appendix I) 
850-mb 00Z 20 August 1965. (Figure 15) 

Four distinct SS ridge lines appear: (1) from just north of Newfoundland 
through New England south to West Virginia and west to Kansas; (2) from 
North Dakota across western Canada to the Pacific; (3) extreme northern 
Canada; and (4) western Alaska. The position of (1) seems well defined by 
a dense data network in its northeastern portion where the position is 
dictated by reports at only five stations: 816, 815, 807, 712, and 600. 
Since these five stations cover 250,000 square miles, location northeast 
of station 712 is questionable. The position at 712 looks very good, since 
a strong maximum occurs there at 700-mb, indicating that the 700-mb front 
is over the station and, therefore, the 850-mb front should be southeast of 
712, as plotted. Along the eastern seaboard the ridge line lies within the 
.150 contour, and its western end is anchored by high readings at stations 
456 and 451. This ridge verifies well with the NMC front, with differences 
of less than fifty miles along 50% of its route. The position of the ridge 
line seems to split the difference between the NMC front and the GG® ridge. 
The ridge over western Canada is based on just a few readings, and would 
not have been drawn in except for the reading at station 867. Itis in- 
teresting to note that the ridge line could have been drawn from station 945 
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north to the strong maximum over station 938. This points out the subjec- 
tivity inherent in hand analysis in sparse data areas. The ridge in northern 
Canada is based entirely on the enormous maximum over station 938. The 
ridge is well defined in the heart of a baroclinic zone. Lack of data 
prevented continuing this ridge farther east. Negative stability layers are 
found at several western stations. These negative values result from one 
significant layer showing % 7 ¥d , a not uncommon occurrence over the 
desert areas in Summer at 850-mb. It is also interesting to note high SS 
along the California coast. This is probably the trade- wind inversion 
(Figure 3) affecting stability at 850-mb. Finally, many readings in the 
Rocky Mountain area are missing because of altitude. 

700-mb 00Z 20 August 1965. (Figure 16) 

The pattern here shows four ridges. The ridge across the eastern 
United States is the same ridge which appeared on the 850-mb chart and 
Shows excellent vertical consistency, in that it trails the 850-mb ridge 
by about 100 miles. If one considers the slope of a cold front to be about 
1/100 between 850 and 700-mb, the ridge position looks reasonable. 

SS trails GG@ over the eastern one third of its path and the distance between 
the two ridges broadens to 150 miles over Missouri and Iowa, but northwest 
of South Dakota the two ridges are in very close proximity. The ridge in 
northeastern Canada is based on one large reading at station 090. This 

ridge could just as well have been omitted, since there are no hints as to 
Orientation and direction. A ridge over the Pacific Northwest is well defined, 
but does not lie in a baroclinic zone, and does not have a counterpart at 


the 850- and 500-mb levels. A ridge across Alaska east and south into 
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central Canada compares fairly well with a GG@ ridge. The SS ridge is 
much longer on its southern end, however, based on high SS over stations 
945 and 848. The vertical consistency of this ridge is not outstanding, 
but this could well be a function of the sparse data area. The western 
portion of this frontal structure which lies in a baroclinic zone, shows 
vertical consistency, but over central Canada there is no baroclinic zone, 
and the front meanders noticeably. A general feature of this chart is the 
flat gradient over the southeast, south central, and Rocky Mountain states. 
SS is not high except over station 304, and since SS increases to the north, 
no ridge was drawn even though the 500-mb analysis shows a ridge in the 
vicinity. 
500-mb 00Z 20 August 1965. (Figure 17) 

The ridge pattern at 500-mb shows four ridges. The 500-mb ridge 
over the eastern United States shows high positive correlation with GG@ 
and with the SS ridge at 850 and 700-mb from Newfoundland to Ohio. 
Northwest of Ohio SS moves steadily away from GG@, taking a more northerly 
route out of the baroclinic zone. It may well be that the SS ridge does follow 
GG8 and maintains vertical continuity with the ridges at lower levels, but 
high readings over stations 747 and 645, plus low readings at 945 and 879, 
make it difficult to draw a ridge across southwestern Canada. A well defined 
ridge over the southern states is not considered frontal, since it is not 
vertically consistent with a lower level front and does not lie in a baroclinic 
zone. The same comments can be made about a ridge from Colorado to 


Washington. However, the ridge over Alaska into Canada is vertically 
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consistent with ridges mentioned at lower levels, and lies in a baroclinic 
zone. Whether this ridge actually ties into the ridge over the midwestern 
United States is conjecture. It is not Surprising that certain difficulties 
are encountered at 500-mb. Some cold fronts may not extend to this 
altitude, or the character of the front may change. As pointed out by 
Gates, the value of static stability may begin to reduce at this level, a 
fact which may affect the pattern. 

300-mb 00Z 20 August 1965. (Figure 18) 

This chart appears to be of little value in location of ridge lines 
which represent fronts. It does seem that SS values are lower at this level, 
as evidenced by the extensive flat gradient over the entire United States, 
with average values of .060. There are three definable ridge lines on the 
chart, but none are associated with lower level ridges and only the one 
over the desert southwest lies in a baroclinic zone. 

It seems logical that fronts, if they exist at all, will be hard to find 
at 300-mb. Since data are sometimes taken up to 250-mb, some of the 
higher stability values may be caused by the tropopause instead of by 
frontal layers. Thus the value of this particular system at 300-mb is 
questionable. For example, analysis of the soundings for stations 906, 
909, and 913 show the tropopause between 250:- and 300-mb. This stable 
zone gives high values of SS and a false indication of a frontal zone. 
850-mb 1200Z 20 August 1965. (Figure 20) 

The major ridge along the eastern seaboard shows time continuity in 
that it has progressed at 15 knots over the eastern states and 25 knots over 


the southern states. The ridge remains very close to the NMC front through 


48 


the area, with a maximum Separation of 100 miles near Cape Hatteras. 

The GG®@ ridge also correlates well, except over Newfoundland and a 
questionable area over the southeastern states. In this area, GG@ turns 
north as the baroclinic zone separates into two branches. The SS ridge 
continues across northern Texas, then northwest into the Pacific, anchored 
by a series of very high readings across the western United States. Another 
secondary ridge, also in a baroclinic zone, extends from Delaware west to 
South Dakota. This matches well with GG@, but has only a small NMC 
counterpart over South Dakota and has no time continuity with any ridge 

on the 00Z chart. This ridge could be the result of a subsidence inversion 
in the cool air behind the front. The ridge over Alaska and northwest 
Canada is essentially stationary from 00Z, and still lies in a baroclinic 
zone. This ridge is more extensive than at 00Z, however. A new ridge over 
extreme northern Canada appears very strong, but there is no baroclinic 
zone in the area. 

An interesting feature of this chart is a large mass of air with low 
stability which dominates central Canada. This mass will be found to 
persist at all levels, and seems to be a feature of the atmosphere behind 
an advancing cold front. 
700-mb 1200Z 20 August 1965. (Figure 21) 

The ridge line along the eastern. seaboard compares very closely with 
GGB6 over its entire path from southeast Canada to Montana, but west of 
Delaware it is closely related to the ridge at 850-mb which was considered 
the secondary ridge. The ridge across the southern states which was strong 


at 850-mb is still strong at 700-mb, but is not continuous across North 


49 


Carolina and Virginia. Considering the high reading at station 317, a 

higher reading at station 402 would have allowed a ridge to intersect the 
strong ridge near station 506, an event which would have given out- 
standing correlation with the SS ridge as shown at 850-mb. The ridge over 
extreme northern Canada has vertical continuity with the 850-mb ridge, but 
still does not lie in a baroclinic zone. The ridge over northwest Canada and 
Alaska looks very realistic over its western portion, but the eastern half is 
questionable because of a lack of data. 

500-mb 1200Z 20 August 1965. (Figure 22) 

The ridge across the northern United States is well defined and has 
consistency with its GG@ counterpart. It also remains imbedded in the 
baroclinic zone. The ridge over the southeast is not in a baroclinic zone, 
and appears to be non-frontal. The ridge over the southwest desert is ina 
baroclinic zone and compares well with GG@ in the area. The ridge over 
Alaska and northwest Canada is still stationary over its western portion, 
but enters a strong north-south baroclinic zone over western Canada, and 
lies very close to the GG@ ridge in the area. Another short ridge over central 
Canada is based on two readings and appears non-frontal. 

Again at 500-mb a large area of low stability exists over eastern 
Canada behind the ridge. 
300-mb 1200Z 20 August 1965. (Figure 23) 

The chart is characterized by small SS over the entire United States 
with the exception of two readings of .15 at stations 583 and 576 which 
define a short ridge over the Nevada, Idaho, Wyoming area. Higher values 


exist over Canada, with a ridge stretching from the Great Lakes north, but 
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the tropopause is responsible for high values of SS at stations 109, 722, 
and 943. There is little data to support placing this ridge entirely inside 
the baroclinic zone outlined by GG@, but contours of SS are generally 
parallel to the ridge and trough lines of GG@. As indicated in the dis- 
cussion of SS at 300-mb for 00Z, there appears little resemblance between 


the pattern at this level and those below. 
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APPENDIX IV 
SS Analysis by Method C 
Comparison of SS with the FNWF Static Stability (0°) Analysis 

Technique and Criteria 

FNWF's hemispheric printout of 0 for 00Z 20 August 1965 displays 
areas of positive and negative o (relative to standard reference values) /6/. 
In some cases these areas are large, with a noticeable gradient, while in 
Other cases small pockets of positive or negative values show up. 

The printouts for the lower four layers were analyzed for ridge lines 
of 0. This analysis was somewhat subjective since no specific guide 
lines were available. Whenever possible, positive areas were connected 
by a ridge line, even though a small area of negative values separates the 
positive areas. When a small positive pocket was completely surrounded 
by negative values, no ridge was drawn. The ridge lines themselves were 
smoothed (i.e., avoiding kinks and sharp corners), whenever feasible, in 
Order to present a more reasonable analysis. No distinction was made 
between a narrow, well defined, positive zone and a wide, poorly defined 
one. No difference was noted between low positive and high positive values 
either, as long as they represented relative maxima. 

These ridge lines were transferred to Figures 15-18 and compared with 
SS ridges. 
850 mb 00Z 20 August 1965. (Figure 15) 

95 at 850-mb is compared with the FNWF O for the layer 1000-775-mb. 
FNWF's O isa medncentc ridge line, which moves in and out of the baro- 


clinic zone. If one assumes that the GG@ ridge and the NMC front are 
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properly located, then O is a rather poor representation of a frontal zone, 
especially over the middle west. The SS ridge, on the other hand, is 
smoother, comparing closely with GG@ and NMC and lying almost entirely 
within the baroclinic zone. SS and O have a close comparison only over 
the middle Atlantic states. One factor tending to produce these variations 
is that O considers the atmospheric structure near the surface -- that is, 
near 1000-mb. Thus, the static stability picture may be confused by layers 
of air near the surface which represent local terrain features. For example, 
the wide swing of O&O out over the Atlantic Ocean may have been caused by 
relatively stable air over the ocean east of Boston -- a not uncommon weather 
feature in summer. On the other hand, SS is computed for the layers around 
850-mb, and this altitude may be above low level inversions. The other 
two oO ridges over the western United States and northern Canada compare 
poorly with GG@ and SS. 
700-mb 00Z 20 August 1965. (Figure 16) 

Considering the ridges over the eastern United States, SS and & 
parallel each other with SS generally 120 miles on the warm air side of & 

The o ridge over central and western Canada is in an area of low 
SS and does not compare well, except over northwestern Canada and Alaska. 
The two small © ridges over the southeast United States and eastern 
Canada are not verified with GG@ or SS ridges. 
500-mb 00Z 20 August 1965. (Figure 17) 

At this level both ridges have less correlation with the GG@ ridge 
and the baroclinic zones than in lower layers, but GO is particularly poor, 
with only about 25% lying in a baroclinic zone and almost no O ridge to 
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compare with GG@. There are at least two locations on the chart where 
ridges cross GG@ ridges at right angles. O does have a slightly better 
positive correlation with SS, possibly because the two quantities have 
some Similarities. 

300-mb 00Z 20 August 1965. (Figure 18) 

Neither O nor SS ridges can justifiably be called good representations 
of frontal phenomena at this level, on the basis of this chart. As stated 
elsehwere in this project, neither ridge seems to give any meaningful 
information at 300-mb. 
850-mb 12002 20 August 1965. (Figure 20) 

The O ridge over the northeastern United States parallels SS and GG8@ 
ridges, generally about 100 miles on the warm side of both. However, over 
the middle western states the O ridge is weak and shows poor comparison 
with extensive GG8@ and SS ridges in the same general area. 

The GO ridges over the far western states and western Canada cross 
the SS ridges at right angles in at least two areas, and there is little basis 
for comparison. The same is true of the ridge over northern Canada when 
compared with GG@ and SS ridges over Alaska. 
700-mb 1200Z 20 August 1965. (Figure 21) 

The highest correlation between SO and SS exists at this level. 

This is a repeat of the situation which existed at 00Z. The oO ridge 
meanders across the SS ridge, while crossing the United States, but 


generally remains within 100 miles of the SS ridge. 
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The O ridge over southwestern Canada is relatively weak and has 
no correlation with a SS ridge except over Alaska. 
500-mb 12002 20 August 1965. (Figure 22) 

At this level O andSS ridges compare well over the northeastern 
United States, but the only other areas where comparison is possible are 
over Alaska, a small area over California, and over the central United 
States. In each case, the O ridge leads the SS ridge, about 100 miles on 
the warm side of the SS ridge. 
300-mb 1200Z 20 August 1965. (Figure 23) 

The major O ridge at this level compares fairly well with the SS ridge, 
and very well with the GG6@ ridge over Canada. 

The analysis over the western states has no SS ridges with which to 
compare the O ridges. Thetwo oO ridges are weak, however, and 
their existence is somewhat subjective. 

Appendix IV pointed out the dominance of low values of SS at 300-mb, 
and the O printout is also dominated by low values, the majority of which 
are negative. Both of these findings agree with the conclusions of Gates, 


who found lower values of static stability above 500-mb VE) 
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